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The current status of the muon anomalous magnetic moment problem is 
briefly presented. The corrections to muon anomaly coming from the effects 
of hadronic vacuum polarization, Z*77* effective vertex and light-by-light 
scattering are estimated within the instanton model of QCD vacuum. 

PACS numbers: 13.40.Em, 14.60.Ef 

1. Muon AMM: experiment vs theory. 

The study of anomalous magnetic moments (AMM) of leptons, a = 
{gs — 2)/2, have played an important role in the development of the standard 
model (SM). At present accuracy the electron AMM due to small electron 
mass is sensitive only to quantum electrodynamic (QED) contributions. The 
theoretical error ^ is dominated by the uncertainty in the input value of 
the QED coupling a = e2/(47r). Thus, the electron AMM provides the best 
observable for determining the fine coupling constant 

= 137.035 998 83(51). (1) 

Compared to the electron, the muon AMM has a relative sensitivity 
to heavier mass scales which is typically proportional to (m^/me)^.^ At 
present level of accuracy, the muon AMM gives an experimental sensitivity 
to virtual W and Z gauge bosons as well as a potential sensitivity to other, 
as yet unobserved, particles in the few hundred GeV/c^ mass range. The 
muon AMM is known to an unprecedented accuracy of order of 1 ppm. The 

* Presented at XLV Cracow School of Theoretical Physics, Zakopane, Poland June 3- 
12, 2005. 

^ The r-lepton AMM due to t's highest mass is the best for searching for manifestation 
of effects beyond SM, however, r-lepton is short living particle, so it is not easy to 
make experiment with good enough accuracy. 
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latest result from the measurements of the Muon {g — 2) collaboration at 
Brookhaven is j2] 

^1(9,-2) = 11 659 208 (6) • W~'\ (2) 

which is the average of the measurements of the AMM for the positively 
and negatively charged muons (Fig. In future, one expects to achieve 
more than a factor of 2 reduction in uncertainty in planning BNL E969 
experiment [3] and even more precise g-2 experiment is discussed in J-PARC 
with the proposal to reach a precision below 0.1 ppm 0]. 
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Fig. 1. Measurements of by E821 (g-2) Collaboration with the SM predictions 
(Theoretical point based on usage of e+e~ ^hadrons annihilation data is 
raised up after recent analysis by SND collaboration 5 ). 

The standard model prediction for consists of quantum electrody- 
namics, weak and hadronic contributions (schematically presented in Fig. 
12]). The QED and weak contributions to have been calculated with great 
accuracy P 

^QED ^ ggg 471.935(0.203) • 10"^° (3) 

and 

a^w ^ 15.4(0.3) . 10-10. (4) 

The uncertainties of the SM value for (Fig. ^ are dominated by the 
uncertainties of the hadronic contributions, a^*'^""^, since their evaluation 
involve quantum chromodynamics (QCD) at long-distances for which per- 
turbation theory cannot be employed. Under assumption that at reached 
scales there are no New Physics effects one may estimate the hadronic part 
of the muon AMM by subtracting the QED and EW contributions from the 
experimental result 



^Strong(Exp) ^ ^21(6) • 10" 



(5) 
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Fig. 2. The standard model contributions to the niuon anomalous magnetic mo- 
ment. 

Below we discuss with some details theoretical status of hadronic contri- 
butions. First, we discuss the phenomenological estimates of the leading of 
order (LO) hadronic corrections based on usage of inclusive e"'"e~ ^hadrons 
and hadronic r decays data. Then, one evaluate the hadronic corrections of 
leading and next-to-leading (NLO) order to muon AMM within the instan- 
ton liquid model of QCD vacuum (ILM). 

2. Phenomenological estimates of the LO hadronic contributions 

to muon AMM 

The LO contribution to the muon AMM comes from the hadronic vac- 
uum polarization (Fig. and NLO corrections consisting of contributions 
which are the iteration of the LO term (Fig. Et) plus the independent con- 
tribution from the light-by- light scattering process (Fig. In absolute 
value the LO and NLO terms differ by one order of magnitude, but the 
theoretical accuracy of their extraction is comparable and dominates the 
overall theoretical error of the SM calculations. All hadronic contributions 
are sensitive to the low energy physics and there are no rigorous theoretical 
methods based on first principles for the calculations. Thus, to confront 
usefully theory with the experiment requires a better determination of the 
hadronic contributions. 

The LO correction to muon AMM, aji^^ , is due to the hadronic photon 
vacuum polarization effect in the internal photon propagator of the one-loop 
diagram (Fig. Using analyticity and unitarity (the optical theorem) 
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o-i^^ can be expressed as the spectral representation integral [71 |H1 

a'r''' = f dt]Kit)p^^\t) , (6) 

which is a convolution of the hadronic spectral function py (t) with the 
known QED kinematical factor 

Jq + [1 - x)t / mf^ 

where is the muon mass. The QED factor is sharply peaked at low t 
and decreases monotonically with increasing t. Thus, the integral defining 
Qhvp(i) sensitive to the details of the spectral function (t) at low 
invariant masses. 
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Fig. 3. Spectral density as measured from e^e -^hadrons annihilation, R{s) — 
At present there is no direct theoretical tools that allow to calculate the 

/TT\ 

spectral function with required accuracy. Fortunately, py (t) is related to 
the total e^e" ^ 7* ^ hadrons cross-section a{t) by (mg — > 0) 

^e+e-^hadrons(^) = iW^p'y'^ (t) , (8) 

and this fact is normally used to get quite accurate estimate of a^^^^\ 
The condensed form accumulating the data of different experiments on the 
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hadronic e^e annihilation is presented in Fig. IHI Moreover, high precision 
inclusive hadronic r decay data 9, 10 , Jj] are used in order to improve the 
determination of aj^^''^^''. This is possible, since the vector current conser- 
vation law relates the 1=1 part of the electromagnetic spectral function 
to the charged current vector spectral function measured va. t ^ v +non- 
strange hadrons. At present, it is found consistence within the experimental 
errors between e^e~ and r data ^5. (see Fig. All these allows to reach 
during the last decade a substantial improvement in the accuracy of the 
contribution from the hadronic vacuum polarization. 
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Fig. 4. The ratio of the e+e^ tt+tt^ cross section calculated from the t~ —* 
TT^TT^Vr decay spectral function measured in [111 1^ to the isovector part of the 
e+e" 'K'^'K~ cross section measured by SND (from ISI). The shaded area shows 
the joint systematic error. 

About 91% of aJi^P^^^ comes from t < (1.8 GeV)^, while 73% of the 
corresponding integral is covered by final 27r state. The most recent es- 
timates of the dispersion integral for the 27r-channel in the energy range 
0.39 < t,r < 0.97 GeV^ which are based on the e~^e~ experimental results 
are following 

aj;" = (378.6 ± 5.0) • 10"^° CMD2 (2003) [12], 
aj;" = (375.6 ± 5.7) • 10"^° KLOE (2004) [13], 
aj;" = (385.6 ± 5.2) • 10"^° SND (2005) [H]. (9) 

The contributions of hadronic vacuum polarization at order quoted in the 
theoretical articles on the subject are given in the Table 1. However, these 
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analysis do not take into account recent SND data which alone may increase 
the estimates based on e^e~ annihilation by approximately (7-^ 10) • lO"^'^ 
(see ©) making e"'"e~ and r data analysis more consistent from one side 
and more close to experimental result from other one. 
Table 1. 

Phenomenological estimates and references for the leading order hadronic 
photon vacuum polarization contribution to the muon anomalous magnetic 
moment based on e^e~ and r data sets. 





e^e il4i 


m 


e+e il^ 


e+e ;16i 


ril6i 




696.3 ±9.8 


711.0 ±8.6 


694.8 ± 8.6 


693.5 ±9.0 


701.8 ±8.9 



The higher order hadronic corrections to are schematically presented 
in Figs. 121 e and g. These diagrams, like leading order contribution, cannot 
be calculated in perturbative QCD, but part of them may be estimated 
with help of experimental data on inclusive hadronic e^e~ annihilation and 
r decays as [T7] ^ 

„hvp{2) ^ _io.l(0.6) • 10-^°. (10) 

This, however, not the case for the so called light-by-light contribution, 
^h. LxL^ (Fig. [2t) where one needs to explore the QCD motivated ap- 
proaches. The latter has been estimated recently using the vector meson 
dominance model supplemented by perturbative QCD constraints [20] 

^h. LxL ^ 13 6(2.5) . 10-1°. (11) 

The agreement between the SM predictions and the present experimen- 
tal values is rather good. There is certain inconsistencies in use of different 
sets of experimental data based on the e^e~ and r processes in evaluations 
of the LO hadronic contribution to the muon AMM. The analysis based 
on the T decay data and recent e+e~ data from SND collaboration |H] pro- 
vide the SM results which are in good agreement with the experimental 
one. The results based on the e^e~ data published by the CMD and 
KLOE 13 data support bigger difference between SM prediction and (g-2) 
Collaboration result. Theoretically, the r decay data is found [U to be 
more compatible with expectations based on high-scale as{Mz) determi- 
nations; the electroproduction data (CMD, KLOE), in contrast, requires 
significantly lower as{Mz)- The results favor determinations of the leading 
order hadronic contribution to which incorporate hadronic r decay data 

^ The second order kernel K^'^\t) has been evaluated in analytical form in |18|. For 
new formulation of the problem of vacuum polarization effects in higher order con- 
tributions to {g — 2) see |19| . 
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over those employing electroproduction data only, and hence suggest a re- 
duced discrepancy between the SM prediction and the current experimental 
value of a^. 



3. The Adler function and a^^'P^^K 

Recently, the isovector vector (V) and axial- vector (A) spectral functions 
have been determined separately with high precision by the ALEPH j^I 
and OPAL jJO] collaborations from the inclusive hadronic r-lepton decays 
(r — > 1/,-+ hadrons) in the interval of invariant masses up to the r mass, 
< s < m^. The vector spectral function measured by ALEPH is shown in 
Fig. 13 It is important to note that the experimental separation of the V 
and A spectral functions allows us to test accurately the saturation of the 
chiral sum rules of Weinberg-type in the measured interval. On the other 
hand, at large s the correlators can be confronted with perturbative QCD 
(pQCD) thanks to sufficiently large value of the r mass. 
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Fig. 5. The isovector vector spectral function (l/47r^)pv(t) from hadronic r- decays 
P|. The dashed line is the asymptotic freedom prediction, (l/47r^). 

Model estimates of the light quark strong sector of the standard model 
will be discussed in the chiral limit, when the masses of it, d, s light quarks 
are set to zero. In this approximation, the V and non-singlet A current- 
current correlation functions in the momentum space (with -g2 = Q2 > Q) 
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are defined as 

Km =ijd^x e'^^-nj^fix) = [{q^q, - g,,q^) Ylj{Q^ , (12) 

n;^,(x) = (0|T{j^(x)J.(0)t}|0), 

where in the local theory the QCD V and A currents for light quarks are 
defined as 

JiM = qiiiVq, J\ = qixi5M, (13) 

the quark field has color (i) and flavor (/) indices, A^^^ = T3 is the isospin 

matrix of the axial current, and V = \{\ + t^) is the charge matrix. The 
momentum-space two-point correlation functions obey (suitably subtracted) 
dispersion relations, 

f 00 r7s 1 

where the imaginary parts of the correlators determine the spectral functions 

pj(s) = 47rlmnj(s + i0). (15) 

Instead of the polarization function it is more convenient to work with the 
Adler function defined as 



n (n^\ n2^nj(Q2) ^ r 



Then, it is possible to express a^^^^^ given by © in terms of the Adler 
function by using the integral representation j22j 

ar(^) = t dx^'-^^^'-^^Dy ( ^ml] , (17) 



^ 3 Jo X \l-x 

where the charge factor YlQf — ^/^j i = u,d, s, is taken into account. The 
bulk of the integral in ()17|) is governed by the low energy behavior of the 
Adler function DyiQ"^). 

The behaviour of the correlators at low and high momenta is constrained 
by QCD. In the regime of large momenta the Adler function is dominated 
by pQCD contribution supplemented by small power corrections 

(18) 
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where the pQCD contribution with three-loop accuracy is given in the chiral 
hmit in MS renormalization scheme by ,23. .2^^ 



nPQCD 




^ I (/^') ^ 



vr 



vr 



+ 



(19) 



+ 
where 



(2F2/3o+/3i)ln^ 



+ ln- 



^2 



vr 



/?0 = 



Pi 



1 



19 

51--n, 



F2= 1.98571 - 0.115295 n/ , ^3= -6.63694 - 1.20013 n/-0.00518n^ , 
with as((5^) being the solution of the equation 



vr 



A 



In 



+ 



In 



A2 ■ 



(20) 



In H18|) along with standard power corrections due to the gluon and quark 
condensates we include the unconventional term suppressed as, ~ 1/Q^. 
Its appearance was augmented in 0U| and also found in the ILM |27j . 
In the 1ow-Q2 limit it is only rigorously known from the theory that 



Dv{Q^^ 0) = Q^D'y{Q)+0{Q^ 



(21) 



It is clear (see also Fig. [71) that the Adler function is very sensitive to 
transition between asymptotically free (almost massless current quarks) re- 
gion described by (|18jl . (|19|) to the hadronic regime with almost constant 
constituent quarks where one has 1)21(1 . 

To extract the Adler function from experimental data supplemented by 
QCD asymptotics ((TH|) . (fT^ we take following an ansatz for the hadronic 
spectral functions in the form 



PJ it) 



ALEPH 
P3 



where 



47r2 



it) 



^PQCD 



12l7r2 fas{t) 



48 



TT 



(22) 



(23) 



and find the value of continuum threshold sq from the global duality interval 
condition: 



so 



dtp-, 



ALEPH 



it) 



so 



dtpf^^ (t) . 



(24) 
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Using the experimental input corresponding to the r-decay data and the 
pQCD expression 



1 /"^o 



N, 



as(so) 



127r2 



[F3 + (2F2/3o + /3i)+2/32 



as(so) 



TT 



+ 



(25) 



as(so) 



vr 



one finds (Fig. that matching between the experimental data and theo- 
retical prediction occurs approximately at scale sq ~ 2.5GeV^. 
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Fig. 6. The integral, Eq. it^ . versus the upper integration limit, sq, for the V 
spectral density. The integral of the experimental data corresponds to solid line 
and the pQCD prediction is given by the dashed line. 

The vector Adler function (jlbj) obtained from matching the low momenta 
experimental data and high momenta pQCD asymptotics by using the spec- 
tral density (|22|) is shown in Fig. Q where we use the pQCD asymptotics 

of the massless vector spectral function to four loops with = 372 

MeV and choose the matching parameter as sq = 

2.5 GeV"^ Admittedly, 

in the Euclidean presentation of the data the detailed resonance structure 
corresponding to the p and oi mesons seen in the Minkowski region (Fig. 

is smoothed out, hence the verification of the theory is not as stringent 
as would be directly in the Minkowski space. 

The phenomenological definition of the Adler function can be used for 
evaluation of the LO contribution to AMM. Below we are going to discuss 
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Fig. 7. The Adler function from the ILM contributions: dynamical quark loop 
(short dashed) , quark + chiral loops + vector mesons (full line) versus the ALEPH 
data (dashed). The dash-dotted line is the prediction of the constituent quark 
model (extended NJL) and the dotted line is the asymptotic freedom prediction, 



the QCD model based definition of the Adler function within the instanton 
liquid model [2^1 • Next two sections we devote to the formulation of the 
gauged instanton liquid model pUj . 



Hadronic corrections to AMM are represented as the convolution inte- 
grals of some known kinematical functions times the amplitudes involving 
low energy quark processes. To study nonperturbative effects of these am- 
plitudes at low momenta one can use the framework of the effective field 
model of QCD. In the low momenta domain the effect of the nonperturbative 
structure of QCD vacuum become dominant. Since invention of the QCD 
sum rule method based on the use of the standard OPE it is common to 
parameterize the nonperturbative properties of the QCD vacuum by using 
infinite towers of the vacuum expectation values of the quark-gluon opera- 
tors. From this point of view the nonlocal properties of the QCD vacuum 
result from the partial resummation of the infinite series of power correc- 
tions, related to vacuum averages of quark-gluon operators with growing 



l/4^2_ 



4. The instanton effective quark model 
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dimension, and may be conventionally described in terms of the nonlocal 
vacuum condensates [311 132j . This reconstruction leads effectively to nonlo- 
cal modifications of the propagators and effective vertices of the quark and 
gluon fields at small momenta. 

The adequate model describing this general picture is the instanton liq- 
uid model of QCD vacuum describing nonperturbative nonlocal interactions 
in terms of the effective action |29| . Spontaneous breaking of the chiral sym- 
metry and dynamical generation of a momentum-dependent quark mass are 
naturally explained within the instanton liquid model. The nonsinglet and 
singlet V and A current-current correlators and the vector Adler function 
have been calculated in |271 1331 134j in the framework of the effective chi- 
ral model with instanton- like nonlocal quark-quark interaction j^. In the 
same model the pion structure function and the pion transition form 
factor normalized by axial anomaly has been considered in |^ for arbitrary 
photon virtualities. The nonperturbative properties of the triangle diagram 
has been thoroughly discussed in P7l 05] . 

We start with the nonlocal chirally invariant action which describes the 
interaction of soft quark fields [HO] 

4 

S = Jd^^ Qlix) [i^D^ - TTif] qi{x) + ^Gp j d^X j ^ d^x„/(x„) 

(26) 

■ \Q{X - xi,X)TpQ{X, X + X3)g(X - X2,X)TpQ{X, X + X4)] , 

where L*^ = — iVf^ (x) — i^5A^ (x) and the spin-flavor structure of the 
nonlocal chirally invariant interaction of soft quarks is given by the matrix 
products^ 

Gp(rp®rp): G(l® l-Fi75T"®i75T"), G' {t'' (g) t"- + i-f5 (g> i-f5) , 

(27) 

where G and G' are the 4-quark couplings in the iso-triplet and iso-singlet 
channels, and r" are the Pauli isospin matrices. For the interaction in the 
form of 't Hooft determinant one has the relation G' = —G. In general due to 
repulsion in the singlet channel the relation G' < G is required. In Eq. ()26() 
Qj = {u, d) denotes the flavor doublet field of dynamically generated quarks. 
The separable nonlocal kernel of the interaction determined in terms of form 
factors f{x) is motivated by instanton model of QCD vacuum. 

In order to make the nonlocal action gauge-invariant with respect to 
external gauge fields V^{x) and A^{x), we define in the delocalized 



^ The explicit calculations below are performed in SUf{2) sector of the model. 
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quark field, Q{x), by using the Schwinger gauge phase factor 



Q{x,y)=Q^x,y)^^, 



(28) 



where P is the operator of ordering along the integration path, with y 
denoting the position of the quark and x being an arbitrary reference point. 
The conserved vector and axial- vector currents have been derived earlier in 

[snuzTiEii- 

The dressed quark propagator, S{p), is defined as 



S-\p) = ip-Mip^), 



(29) 



with the momentum-dependent quark mass found as the solution of the gap 
equation 



M{p')=mf + AGpNfNj\p') J -^J 



(30) 



(31) 



The formal solution is expressed as 

M{p^)=mf + {Mg-mf)f{p^), 

with constant Mq = M(0) determined dynamically from Eq. (|30j) and 
the momentum dependent f{p) is the normalized four-dimensional Fourier 
transform of /(x) given in the coordinate representation. 

The nonlocal function f{p) describes the momentum distribution of 
quarks in the nonperturbative vacuum. Given nonlocality f{p) the light 
quark condensate in the chiral limit, M{p) = Mgf'^{p), is expressed as 



(0 \qq\ 0) 



d^p M{p'^) 



(32) 



Its n-moment is proportional to the vacuum expectation value of the quark 
condensate with the covariant with respect to gluon field derivative squared 
to the nth power 



{0\qD^^q\0) = -N,l 



d'^p 



.2ra 



M(p2 



47r4^ p'^+M^ip'^)' 



(33) 



The nth moment of the quark condensate appears as a coefficient of Taylor 
expansion of the nonlocal quark condensate defined as |SJ 



C{x) 



q (0) P exp 



i / A^{z)dz^ 



q{x] 







(34) 
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with gluon Schwinger phase factor inserted for gauge invariance and the 
integral is over the straight Une path. Smoothness of C{x) near = 
leads to existence of the quark condensate moments in the l.h.s. of for 
any n. In order to make the integral in the r.h.s. of ()33() convergent the 
nonlocal function f{p) for large arguments must decrease faster than any 
inverse power of p^, e.g., like some exponential 

/ (p) ~ exp (— const -p") , a > as —>■ oo. (35) 

Note, that the operators entering the matrix elements in and ()34() 
are constructed from the QCD quark and gluon fields. The r.h.s. of (|33() 
is the value of the matrix elements of QCD defined operators calculated 
within the effective instanton model with dynamical quark fields. Within 
the instanton model the zero mode function f{p) depends on the gauge. 
It is implied [HH ESI that the r.h.s. of corresponds to calculations 
in the axial gauge for the quark effective field. It is selected among other 
gauges because in this gauge the covariant derivatives become ordinary ones: 
D — > (9, and the exponential in (|34j) with straight line path is reduced to unit. 
In particular it means that one uses the quark zero modes in the instanton 
field given in the axial gauge when define the gauge dependent dynamical 
quark mass. The axial gauge at large momenta has exponentially decreasing 
behavior and all moments of the quark condensate exist. In principle, to 
calculate the gauge invariant matrix element corresponding to the of l.h.s. of 
H33|) it is possible to use the expression for the dynamical mass given in any 
gauge, but in that case the factor p'^^ will be modified by more complicated 
weight function providing invariance of the answer^. 

Furthermore, the large distance asymptotics of the instanton solution is 
also modified by screening effects due to interaction of instanton field with 
surrounding physical vacuum [321 140j . To take into account these effects and 
make numerics simpler we shell use for the nonlocal function the Gaussian 
form 

/(p)=exp(V/A'), (36) 

where the parameter A characterizes the nonlocality size of gluon vacuum 
fluctuations and it is proportional to the inverse average size of instanton 
in the QCD vacuum. 

The important property of the dynamical mass H3Ujl is that at low virtu- 
alities its value is close to the constituent mass, while at large virtualities it 
goes to the current mass value. As we will see below this property is crucial 

* If one would naively use the dynamical quark mass corresponding to popular singular 
gauge then one finds the problem with convergence of the integrals in 13311 . because 
in this gauge there is only powerlike asymptotics of AI (p) ~ p^" at large p^. 
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Fig. 8. Diagrammatic representation of the Bethe-Salpeter equation for the quark- 
quark scattering matrix, T, with nonlocal instanton kernel, /. 

in obtaining the anomaly at large momentum transfer. The instanton liq- 
uid model can be viewed as an approximation of large- A'c QCD where the 
only new interaction terms, retained after integration of the high frequency 
modes of the quark and gluon fields down to a nonlocality scale A at which 
spontaneous chiral symmetry breaking occurs, are those which can be cast 
in the form of four-fermion operators H26() . The parameters of the model are 
then the nonlocality scale A and the four-fermion coupling constant Gp. 




5. Conserved vector and axial-vector currents 

The quark-antiquark scattering matrix (Fig. |HJ in pseudoscalar channel 
is found from the Bethe-Salpeter equation as 

with the polarization operator being 

•^^^(9') = / 7?T4 /' (^) /' + 9) Tr [S{k)^,S (k + q) 75] . (38) 
J (27r) 

The position of pion state is determined as the pole of the scattering 
matrix 

det(l - GpJpp{q^))\^,^_^, = 0. (39) 

The quark-pion vertex found from the residue of the scattering matrix is 
{k' = k + q) 

Tl{k,k')=g^,,i^^f{k)f{k')r^ (40) 



with the quark-pion coupling found from 

(41) 



9nq 



dJpp iq^) 
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(a) 



\ 



Fig. 9. Diagrammatic representation of the bare (a) and full (b) quark-current 
vertices. Diagram (c) shows separation of local (fat dot) and nonlocal parts of the 
full vertex. 

where m^r is physical mass of the vr- meson. The quark-pion coupling, gj^q, 
and the pion decay constant, f-^, are connected by the Goldberger-Treiman 
relation, g.,^ = Mq/fj^, which is verified to be valid in the nonlocal model 
j39j . as requested by the chiral symmetry. 

The vector vertex following from the model (|26|) is (Fig. ISH) 



mass (see below (|57() ). q is the momentum corresponding to the current, 
and k [k') is the incoming (outgoing) momentum of the quark, k' = k + q. 

The full axial vertex corresponding to the conserved axial-vector current 
is obtained after resummation of quark-loop chain that results in appearance 
of term proportional to the pion propagator [SUj (Fig. Eb) 




Tl{k,k') = ^^-i^ + 2^^^f{k)f{k') JAPiO) 



1 - GpJpp ((/2) 



+ {k + k )^Jap (0) _ fc2 



(43) 



where we have introduced the notations 




(44) 
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JAp{q^) = 4iVciV/ I ^^^/M{fTq)M{r). (45) 
The axial- vector vertex has a pole at 

2 2 V , „2 

q = -m^ = mciqq) /U 

where the Goldberger-Treiman relation and definition of the quark conden- 
sate have been used. The pole is related to the denominator 1 — GpJpp (g^) 
in Eq. (|43j) . while in denominator is compensated by zero from square 
brackets in the limit — > 0. This compensation follows from expansion of 
J{q'^) functions near zero 

Jpp(g2) = Gp' + m, (qq) - q^g-^ + O {q') , (46) 

JAp{q^ = 0) = M„ Jp{q^ = 0) = {qq) M'K (47) 

In the chiral limit = the second structure in square brackets in Eq. 
disappears and the pole moves to zero. 
Within the chiral quark model |3Uj based on the non-local structure of 
instanton vacuum j^lj the full singlet axial- vector vertex including local and 
nonlocal pieces is given by (in chiral limit) (17 \ 

TfHk.g, k'^k + q) = 1^1, + 75(fc + t');.A^, '-'''*2l(i*"' + (48) 

The singlet current (|48|) does not contain massless pole due to presence 
of the Ua (1) anomaly. Indeed, as — > there is compensation between 
denominator and numerator in H48|l 

l-GJppjq') fl 2 n MQ^ 

z^2 ^^mI ^ ^ 

where /tt is the pion weak decay constant. In cancellation of the massless 
pole the gap equation is used. Instead, the singlet current develops a pole 
at the r]'— meson mass^ 

1 - G'Jpp{q^ = -ml,) = 0, (50) 



^ See footnote 3. Also we neglect the effect of the ajdal-pseudoscalar mixing with the 
longitudinal component of the flavor singlet /i meson. 
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thus solving the Ua^^) problem. Let us also remind that in the instanton 
chiral quark model the connection between the soft gluon and effective quark 
degrees of freedom is fixed by the gap equation. In particular, it means that 
the four-quark couplings G{G') are proportional to the gluon condensate. 

The parameters of the model are fixed in a way typical for effective low- 
energy quark models. One usually fits the pion decay constant, 7,^, to its 
experimental value, which in the chiral limit reduces to 86 MeV J41'. In the 
instanton model the constant, is expressed as 

oo 

„2 f , M'^{u)-uM{u)M'{u)+v?M'{uf 



f^ = ^2jduu , (51) 







where here and below u = k'^, primes mean derivatives with respect to u: 
M'{u) = dM{u)/du, etc., and 

D{k^) =k^ + M^{k). (52) 

One gets the values of the model parameters jSl] 

Mg = 0.24 GeV, Ap = 1.11 GeV, Gp = 27 A GeY-^. (53) 

The coupling G' is fixed by fitting the meson spectrum. Approximately one 
has G' « 0.1 G |Sn]. 

6. Adier function within ILM. 

Our goal is to obtain the vector current-current correlator and corre- 
sponding Adler function by using the effective instanton-like model ()26(1 
and then to estimate the leading order hadron vacuum polarization correc- 
tion to muon anomalous magnetic moment a^. In ILM in the chiral limit 
the (axial-)vector correlators have transverse character |27j 

K iQ') = {a,. - ^) nr iQ') , (54) 

where the polarization functions are given by the sum of the dynamical 
quark loop, the intermediate (axial-)vector mesons and the higher order 
mesonic loops contributions (see Fig. I10|l 

uf"^ (q2) = n?^°°p (q2) + n7-- {q^) + (g^) . (55) 

The spectral representation of the polarization function consists of zero 
width (axial-)vector resonances (n™'^^""'^ {Q"^)) and two-meson states fllj'^""^ 
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H 



■k;K 



Fig. 10. Schematic representation of the vector polarization function (|55|l . 



The dynamical quark loop under condition of analytical confinement has no 
singularities in physical space of momenta. 

The dominant contribution to the vector current correlator at space-like 
momentum transfer is given by the dynamical quark loop which was found 
in \n\ with the result^ 



n 



QLoop 
V 



3 ^ 



(2t:Y D+D^ 



k^k^ 



(56) 



(^M(i) {k+,k^)y {k+k^ - M+M^) - (M^ k^)) 



(1) 



where the notations 



8A^c f d'^k M{k) 
i27rfD{k) 



I 



M' (k) - -A:iM(2) {k, k + Q,k) 



k± = k±Q/2, 



ki 



k^k- 



(k+q) (k^q) 



M± = M(fc±), D±=D{k±), 

are used. We also introduce the finite-difference derivatives defined for an 
arbitrary function F (k) as 



F'^^\k,k') 



F{k') - F{k) 



(57) 

In ()56p the first integral represents the contribution of the dispersive 
diagrams and the second integral corresponds to the contact diagrams (see 
Fig. ^2 and ref. for details). The expression for ng^°°P (Q2) is formally 
divergent and needs proper regularization and renormalization procedures 
which are symbolically noted by [--l^g^ for the divergent term. At the same 
time the corresponding Adler function is well defined and finite. 

Also we have checked that there is no pole in the vector correlator as 
0, which simply means that photon remains massless with inclusion 



Within the context of ILM, the integrals over the momentum are calculated by trans- 
forming the integration variables into the Euclidean space, 
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Fig. 11. The dynamical quark-loop contribution is the sum of dispersive and contact 
terms. In the dispersive diagram F is the bare vertex and F is the total one. 

of strong interaction. In the limiting cases the Adler function derived from 
Eq. in accordance with the first equality of Eq. HI6() satisfies general 
requirements of QCD (see leading terms in (fTH|) . (fT^ . and 

iQ' - 0) = o (q2) , (Q' ^ = + ^ + ^ iQ^") ■ 

(58) 

The leading high asymptotics comes from the ^k^k- — f^^lren ^^'^^ ™ 
H56|l . while the subleading asymptotics is driven by "tachionic" term with 
coefficient [2Z] 

It is possible to integrate Eq. (jSU)) in the dilute liquid approximation, u >> 

OX ^ « 4.7 • 10-=^ GeV^ (60) 

which is close to exact result |27j and phenomenological estimate from |26j . 

In the extended by vector interaction model (|26() one gets the corrections 
due to the inclusion of p and u; mesons which appear as a result of quark- 
antiquark rescattering in these channels 

where By (9^) is the vector meson contribution to quark-photon transition 
form factor 

f d^k f^f^ 

Bv (Q') = 8N,i j ^^^^4 [M+M. - k+k.+ (62) 



+ '-kl (l-M^W {k+,k. 
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and Jy (g^) is the vector meson polarization function defined in (jHEI) witli 
^Ji — {dtJ-f ~ (liilu/o^) lu- As ^ consequence of the Ward-Takahashi identity 
one has By (0) = as it should be. 

To estimate the tt+tt" and K~^K~ vacuum polarization insertions (chiral 
loops corrections) one may use the effective meson vertices generated by the 
Lagrangian 

' " " ^ (63) 



-ie I TT 



By using the spectral density calculated from this interaction: 

4^2 X 3/2 



1 

12 



1 



t 



e{t - Ami) + {t^^K 



one finds the contribution to the Adler function as 



where 



a{t) 



nXLoop 



3 + t 



487r2 



a(^\+a 



t + l 



arctanh 



l + 2t 



2v/t(t+l) 



.vr 



(64) 



(65) 



(66) 



The estimate H65() of the chiral loop corrections corresponds to the point- 
like mesons which becomes unreliable at large t, where the meson form 
factors has to be taken into account. This contribution corresponds to the 
lowest order, O(p^), calculations in chiral perturbation theory (^PT), is 
non- leading in the formal l/A'^c-expansion and provides numerically small 
addition. The higher-loop effects become important at higher momenta. 

The resulting Adler function in ILM is given by the sum of above con- 
tributions 

Dv = (Q2) + D^'^°-^ (Q2) + I)XLoop ^^2) _ 

By using set of parameters found in ILM, , the Adler function in the 
vector channel (|67|) is presented in Fig. Q ^-nd the model estimate for the 
hadronic vacuum polarization to given by (|17|) is 



a 



hvp (1);ILM 



623 (40) • 10 



'10 



(68) 



where the various contributions to a 



hvp (1);ILM 



are 



hvp (l);Qloop 



533-10 



'10 



^hvp (l);Vmesons 13.10 



10 



,hvp (l);xLoop 



77- 



(69) 
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and the error in (|(i8j) is due to incomplete knowledge of the higher order 
effects in nonchiral corrections. One may conclude, that the agreement of 
the instanton model estimate with the phenomenological determinations in 
Table 1 is rather good, but model approach unlikely reaches the required 
by experiment accuracy. Nevertheless, for the higher order hadronic cor- 
rections we are able essentially reduce the theoretical error by using rather 
sophisticated effective quark models. The realistic model calculations are 
a crucial issue in consideration of the NLO hadronic contributions. Repro- 
ducing the phenomenological determination of a^^ it becomes possible 
to make reliable estimates of Aa^^ and a^; LxL_ 

With the same model parameters one also gets the estimate for the 
hadronic contribution to the r-lepton anomalous magnetic moments 

^hvp (1);ILM ^ 3 (0 2) • 10^*^ , (70) 
which is in agreement with phenomenological determination 
hvp(i);exp_/ 3.383(0.111) -10-6, 



' ^ 3.536(0.038) • 10-^ IS], 

and prediction of the gauged nonlocal quark model j44j 
^hvp {i);GNQM^ 3 2(0.1). 10-6. 

Thus, we conclude that the LO hadronic corrections obtained within the 
ILM are in reasonable agreement with the latest precise phenomenological 
numbers. Next, we are going to use the ILM in order to estimate a subset of 
hadronic contributions to the muon anomalous magnetic moment, 



7. VAV correlator and NLO corrections to 

Since discovery of anomalous properties [13 IS] of the triangle diagram 
with incoming two vector and one axial-vector currents |17] many new in- 
teresting results have been gained. Recently the interest to triangle diagram 
has been renewed due to the problem of accurate calculation of higher order 
hadronic contributions to muon anomalous magnetic moment via the light- 
by-light scattering process (Fig. I12j) ''. a^' ^xL^ that cannot be expressed 
as a convolution of experimentally accessible observables and need to be 
estimated from theory. 

The light-by-light scattering amplitude with one photon real and an- 
other photon has the momentum much smaller than the other two, can be 



^ See, e.g., |2UI 14911^0) and references therein. 
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q^=o 

Fig. 12. OPE presentation of the light-by-light scattering as the triangle amplitude 
times the coefficient function 




Fig. 13. Effective Z*^^* coupling induced by a fermion triangle contributing to 



analyzed using operator product expansion (OPE). In this special kinemat- 
ics the amplitude is factorized into the amplitude depending on the largest 
photon momenta and the triangle amplitude involving the axial current A 
and two electromagnetic currents (one soft V and one virtual V). The very 
similar kinematics for the triangle amplitude with quark and lepton inter- 
nal lines also defines a subset of the two-loop contributions to aj^^ via the 
Z*77* effective coupling (Fig. 

The corresponding triangle amplitude, which can be viewed as a mixing 
between the axial and vector currents in the external electromagnetic field, 
were considered recently in [101 EH EHl IIHI • This amplitude can be written 
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as a correlator of the axial current j| and two vector currents jy and 

T^.X = - I dWy e*^-'-'^^ (0| T^Cx) ji(0)}|0) , (71) 



where the currents are defined in (|13j) . with the tilted current being for the 
soft momentum photon vertex. In the specific kinematics when one photon 
(92 = 9) is virtual and another one {qi) represents the external electromag- 
netic field and can be regarded as a real photon with the vanishingly small 
momentum qi depends only on two invariant functions, longitudinal wl and 
transversal wt with respect to axial current index j51j . 



-WL (q^) q^qiq''£:pfiau+ 



+WT {q"^) {q^qiBp^ux - q" qiq2^pticj\ + q^qiq2^ppau^ ■ (72) 

Both structures are transversal with respect to vector current, q^T^^x = 0. 
As for the axial current, the first structure is transversal with respect to 
while the second is longitudinal and thus anomalous. 

In the local theory the one-loop result for the invariant functions wt 
and wl is^ 



i-ioop o i-ioop f da a(l - a) 

='"^ =—L a(l-a)<Z^ + m- ^''^ 



where the factor Nc/2> is due to color number and electric charge. In the 
chiral limit, mj = 0, one gets the result for space-like momenta q (jf' > O) 

WL {q^) = 2WT {q^) = ^. (74) 

The appearance of the longitudinal structure is the consequence of the 
axial Adler-Bell- Jackiw anomaly |451 146j . For the nonsinglet axial current 
A^^^ there are no perturbative [HH and nonperturbative |55| corrections to 

(3) 

the axial anomaly and, as consequence, the invariant function w\^ remains 
intact when interaction with gluons is taken into account. Recently, it was 
shown that the relation 

WLT (g') = WL {q^) - 2wT {q^) = 0, (75) 



Here and below the small effects of isospin violation is neglected, considering m/ 
ruu — md- 



dorokhov printed on February 2, 2008 



25 



which holds in the chiral hmit at the one- loop level (|74j) , gets no perturbative 
corrections from gluon exchanges in the iso-singlet case Nonpertur- 
bative nonrenormalization of the nonsinglet longitudinal part follows from 
the 't Hooft consistency condition i.e. the exact quark-hadron duality 
realized as a correspondence between the infrared singularity of the quark 
triangle and the massless pion pole in terms of hadrons. OPE analysis in- 
dicates that at large q the leading nonperturbative power corrections to wt 
can only appear starting with terms ~ containing the matrix elements 
of the operators of dimension six [SH]. Thus, the transversal part of the 
triangle with a soft momentum in one of the vector currents has no pertur- 
bative corrections nevertheless it is modified nonperturbatively. However, 
for the singlet axial current A^'^'> due to the gluonic Ua (1) anomaly there 
is no massless state even in the chiral limit. Instead, the massive r/' me- 
son appears. So, one expects nonperturbative renormalization of the singlet 
anomalous amplitude at momenta below r]' mass. Below we demon- 

(3) 

strate how the anomalous structure is saturated within the instanton 
liquid model. We also calculate the transversal invariant function wt at ar- 
bitrary space-like q and show that within the instanton model in the chiral 
limit at large all allowed by OPE power corrections to wt cancel each 
other and only exponentially suppressed corrections remain ^37^ i38j . The 
nonperturbative corrections to wx at large q^ have exponentially decreasing 
behavior related to the short distance properties of the instanton nonlocality 
in the QCD vacuum. 

The contribution of Z*77* vertex to the muon AMM a^j^ in the unitary 
gauge, where the Z propagator is i {—g^u + QnQi//fn'z) I {'f ~ "^z)' '^^'^ 
written in terms of wl^t (9^) as 



1 

rriy 

+ 2 2 ^^ 

— 



-(1 + (uj _ \ ^ 

(76) 



where p is the four-momentum of the external muon, = 1.16637(1) • 
10^^ GeV~^ is the Fermi constant obtained from the muon lifetime, mz = 
91.1875(21) GeV, a = 1/137.036 and for the electron neglecting its mass 
one has 

2 

In perturbative QCD with massless quarks the result for the first generation 



WL [e] = 2wT [e] = -^. (77) 



This relation for massive quarks is proved to be valid up to two-loop level |54| . 



26 



DOROKHOV PRINTED ON FEBRUARY 2, 2008 





Fig. 14. Diagrammatic representation of the triangle diagram in the instanton 
model with dressed quark lines and full quark-current vertices (a) ; and part of the 
diagram when all vertices are local one (b). 



[e, M, d] contribution is 



Aa 



EW 



e, M, d] = 0, 



(78) 



due to anomaly cancellation. 



8. V AV correlator within the instanton hquid model 

Our goal is to obtain the nondiagonal correlator of vector current and 
nonsinglet axial- vector current in the external electromagnetic field (Fviy) 
by using the effective instanton-like model H26|) . In this model the VAV 
correlator is defined by (Fig. [rik ) 



T^iyx{qi,q2) = -2Nc 



Tr [T^,{k + qi,k)S{k + qi) 



■Tl {k + q,,k-q2)S{k- q2) {k, k - q^) S {k)] 



(79) 



where the quark propagator, the vector and the axial-vector vertices are 
given by (PUjl , (|^^ and , respectively. The structure of the vector ver- 
tices guarantees that the amplitude is transversal with respect to vector 
indices 

Tf,ux{qi,q2)qi = Tf,^x{qi,q2)q2 = 

and the Lorentz structure of the amplitude is given by H72|) . 

It is convenient to express Eq. H79|) as a sum of the contribution where all 
vertices are local (Fig. 114b ). and the rest contribution containing nonlocal 
parts of the vertices (Fig. ITlk). Further results in this section will concern 
the chiral limit. 
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The contributions of diagram ll4b to the invariant functions at space-hke 
momentum transfer, = are given by 



w 



(loc) 



9g2 



1 



7r2 DlD^ 



A;2_4i^ + 3(M 



w 



(loc) 
LT 



(80) 



(81) 



where we also consider the combination of invariant functions wlt, <|75|1 . 
which show up nonperturbative dynamics. The notations used here and 
below are 



k. 



k, 



k. 



k 



ki 



k+k- 



(k+q) (k-q) 



D± = Dik"^), M± = M{kl] 
At large one has an expansion 



w 



(loc) 



oo 



) 



1 



+ O {q-') 



(82) 



(83) 



It is clear that the contribution (jHUj) saturate the anomaly at large g^. The 
reason is that the leading asymptotics of ()8U() is given by the configuration 
where the large momentum is passing through all quark lines. Then the 
dynamical quark mass M{k) reduces to zero and the asymptotic limit of 
triangle diagram with dynamical quarks and local vertices coincides with 
the standard triangle amplitude with massless quarks and, thus, it is inde- 
pendent of the model. 

The contribution to the form factors when the nonlocal parts of the 
vector and axial-vector vertices are taken into account is given by 



w 



(nonloc) 



ANc f d'^k 1 



-M'^^^\k+,k. 



Skq) 



2 



(84) 



Summing analytically the local H80() and nonlocal (|51)) parts provides us 
with the result required by the axial anomaly j37j 



(85) 
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Fig. 15. Normalized wl invariant function constrained by ABJ anomaly from tri- 
angle diagram Fig. 114b (solid line) and different contributions to it: from local 
part, Fig. 114b . (dashed line), and from the nonlocal part (dash-dotted line). 




Fig. El illustrates how different contributions saturate the anomaly. Note, 
that at zero virtuality the saturation of anomaly follows from anomalous 
diagram of pion decay in two photons. This part is due to the triangle 
diagram involving nonlocal part of the axial vertex and local parts of the 
photon vertices. The result (|85jl is in agreement with the statement about 
absence of nonperturbative corrections to longitudinal invariant function 
following from the 't Hooft duality arguments. 

For WLTi'f') a number of cancellations takes place and the final result 
is quite simple fST 



3q^ 



-Mi 



^ki 



y^M^'^\k+,k-) -2{kq)M'^ 



(1) , 
{k+,k^ 



(86) 



The behavior of WLxiq'^) is presented in Fig. 1171 In the above expres- 
sion the integrand is proportional to the product of nonlocal form factors 
/ / (^^) depending on quark momenta passing through different quark 
lines. Then, it becomes evident that the large q^ asymptotics of the integral 
is governed by the asymptotics of the nonlocal form factor / ((7^) which 
is exponentially suppressed (|35|) . Thus, within the instanton model the 
distinction between longitudinal and transversal parts is exponentially sup- 
pressed at large q^ and all allowed by OPE power corrections are canceled 
each other. The instanton liquid model indicates that it may be possible 
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that due to the anomaly the relation (|75|) is violated at large only expo- 
nentially. 

The calculations of the singlet VAV correlator results in the following 
modification of the nonsinglet amplitudes |38j 



where 



w 



(87) 



9g2 1 - G'Jpp (g2) J ^~Dpr~ 



1 - G'/G 



(89) 



4 {kqY 
3^ 



(kq) 



1,0 

d'WJQ) 

0,8 - / 
0,6 - ,' 
0,4 - i' 
0,2 - 'i 

0,0^ • 1 • 1 • 1 • 1 

12 3 4 

Q(GeV) 

Fig. 16. Normalized Wh invariant function in 
the nonsinglet case (solid line) and singlet case 
(dashed line). 




Q(GeV) 



Fig. 17. Normalized wlt invariant function 
versus Q predicted by the instanton model in 
the nonsinglet case (solid line) and singlet case 
(dashed line). 



Fig. El illustrates how the singlet longitudinal amplitude is renor- 
malized at low momenta by the presence of the U a (1) anomaly. The behav- 
ior of w^^^{q'^) is presented in Fig. El Precise form and even sign of w^^)p{q^) 
strongly depend on the ratio of couplings G' / G and has to be defined in the 
calculations with more realistic choice of model parameters. 
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By using (|7(i|) one finds numerically the result for the first generation 
[e, M, d] contribution 

Aaf^[e,u,d]=-lA8-10-^\ (90) 

which has to be compared with recent numbers —2.02 • 10~^^ [HI obtained 
from simple vector dominance model and —4 • 10~^^ |HZ] calculated in the 
naive constituent quark model. 

The preliminary estimate of the hadronic light-by-light scattering con- 
tribution within the instanton liquid model is 

a^^^x^ = 10.6(1.0) -W-^^, (91) 

which has to be compared with in 13.6(2.5) • 10"^'^ where the simple 
vector meson dominance model has been used. 

9. Conclusions 

We briefly discussed the current status of experimental and theoretical 
results on the muon anomalous magnetic moment. The biggest theoretical 
error is due to hadronic part of AMM. The phenomenological and model 
approaches considered for estimates of leading and next-to-leading order 
hadronic corrections to muon AMM. For the model estimates one has used 
the instanton liquid model of QCD vacuum. We calculated the vector Adler 
function and the nondiagonal correlator of the vector and axial-vector cur- 
rents in the background of a soft vector field for arbitrary space-like mo- 
menta transfer and found the corrections to muon anomaly coming from the 
effects of hadronic vacuum polarization, Z*77* effective vertex and light- 
by-light scattering. 
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